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Abstract. Inspired by Raychaudhuri's work, and using the equation named after him as a basic 
ingredient, a new singularity theorem is proved. Open non-rotating everywhere expanding universes 
with non-vanishing spatial average of the matter variables are severely geodesically incomplete to 
the past. Another way of stating the result is that, under the same conditions, any singularity-free 
model must have a vanishing spatial average of the energy density (and other physical variables). 
This is very satisfactory and provides a clear decisive difference between singular and non-singular 
cosmologies. 
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1. Introduction 

In this paper I would like to present a result which confirms — at least partially — a 10- 
year-old conjecture: reasonable non-rotating cosmological models can be non-singular 
only if they are open and have vanishing spatial averages of the matter and other physi- 
cal variables. This conjecture arose as a result of the interactions and discussions between 
Professor A K Raychaudhuri (AKR from now on) and myself concerning the question of 
the feasibility of singularity-free cosmological solutions, see [35,48]. This is a subject to 
which, as is well-known, AKR made fundamental pioneering contributions [33,34], and I 
influenced in a much more modest and lateral manner [43,46,47]. 

The next section contains a more or less detailed historical review of the antecedents, 
birth, and hazardous life of the conjecture. However, any reader, not interested in this 
historical appraisal, who only wishes to see and learn about the results and the new theorem 
— which are interesting on their own — may go directly to Sects. 3. and 4.. The former 
contains a brief derivation of the celebrated Raychaudhuri equation and some discussion on 
the definition of spatial averages; the latter presents the new theorem and its implications. 
A brief discussion and an Appendix have been placed at the end of the paper. 



*In memoriam Amal Kumar Raychaudhuri (1924-2005) 
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2. History of a conjecture 

This is divided into three parts: the appearance of the first non-singular solution and subse- 
quent developments are summarized in subsection 2.2.1; the importance of the fresh new 
ideas brought in by AKR, and the formulation of the conjecture, are analyzed in subsection 
2. 2.2; finally, the reaction by AKR about the conjecture, the sometimes chaotic discus- 
sions in publications-comments-replies, and some recent developments on the subject are 
reviewed in subsection 2. 2.3. 



2.1 Anadi Vishva 

In 1990 1 published [43] what is considered to be the first "interesting" singularity-free cos- 
mological model. It is a spatially inhomogeneous solution of Einstein's field equations for 
a perfect fluid with a realistic equation of state p = g/3. This solution came as a big sur- 
prise: it was widely believed that a model such as the one in [43] should be singular due to 
the powerful singularity theorems developed by Penrose, Hawking and others [31,22,21]. 
Therefore, the solution caused some impact, see e.g. [26], and some discussions within the 
relativity community [23,44]. 

Shortly afterwards, the geodesic completeness of this solution was explicitly proven in 
[2], along with a list of its main properties. The solution turns out to be cylindrically sym- 
metric, globally hyperbolic and everywhere expanding for half the history of the model. 
It satisfies the stronger energy requirements (the dominant and strict strong energy condi- 
tions) and is singularity-free (see the Appendix for a short summary). A detailed analysis 
of how the model fits in with the general conclusions of the singularity theorems was also 
performed in [2], and the solution was proven to be in full accordance with the theorems: 
in all versions of the theorems at least one of their hypotheses was not satisfied, usually the 
so-called boundary or initial condition, see [47]. 

Almost simultaneously, another paper [41] extended the particular model in [43] to a 
general family of perfect- fluid solutions with an Abelian Gi -symmetry acting on spatial 
surfaces. This family was obtained under the assumption of separation of variables. It con- 
tains a diversity of models, with different properties, in particular a 2-parameter subfamily 
of geodesically complete, singularity-free, solutions — which include the original solution 
in [43] — as well as other relevant solutions such as the one in [10], see also [51]. There 
were also solutions with timelike singularities, and solutions with no singularities in the 
matter variables but with timelike singularities in the Weyl conformal part of the curvature 
[41]. All such behaviour was somehow unexpected. As a matter of fact, one could prove 
that, within the subfamily of solutions in [41] with a barotropic equation of state and no 
initial singularity, the singularity-free solutions lie precisely at the boundary separating the 
solutions with only timelike Weyl singularities from the solutions with timelike singular- 
ities in both Weyl and matter variables. This was not very encouraging, since one could 
suspect the instability and the zero measure of the non-singular models. 

Several generalizations of these solutions were later found by allowing the fluid to have 
heat flow [30] and scalar fields [7]. It was claimed in [9,7] that a particular "inhomogenei- 
sation" procedure (by requiring the separation of variables) of the Robertson- Walker open 
cosmological models would lead to the non-singular models of [43,41]. This was later pur- 
sued in a review published by Dadhich in [3], where he intended to prove the uniqueness 
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of those non-singular models. 

Up to this point, all known non-singular solutions were (i) diagonal, that is to say, 
with the existence of a global coordinate chart adapted to the perfect fluid (so-called "co- 
moving") such that the metric takes a diagonal form, (ii) cylindrically symmetric, and (iii) 
separable in comoving coordinates. The last feature meant that the metric components 
could be written as the product of a function of the separable time coordinate times a func- 
tion of the separable radial coordinate. The first example of a non-diagonal non-singular 
perfect fluid model was presented in [27], though the solution had been previously pub- 
lished in quite a different context in [25]. This was a solution for a cylindrically symmetric 
stiff fluid (the equation of state is p = q) or equivalently, for a massless scalar field. It was 
also separable but could be almost immediately generalized to a family of non-separable, 
non-diagonal, non-singular stiff fluid solutions in [20]. Other non-singular solutions fol- 
lowed, see e.g. [28,29,47]. 

It was then shown in [46] that the general family in [41] as well as the whole class 
of Robertson- Walker cosmologies belong to a single unified wider class of cylindrically 
symmetric 1 , separable, diagonal (non-necessarily perfect) fluid solutions. Those depended 
on one arbitrary function of time — essentially the scale factor — and four free parameters 
selecting the openness or closeness of the models, the anisotropy of the fluid pressures, 
or the anisotropy and spatial inhomogeneity of the models. The physical properties of 
this general class were analyzed in detail [46,47], in particular the deceleration parameter, 
leading to natural inflationary models (without violating the strong energy condition), and 
the generalized Hubble law. The possibility of constructing realistic cosmological models 
by "adiabatically" changing the parameters in order to start with a singularity-free model 
which at later times becomes a Robertson- Walker model was also considered in [46] and 
in section 7.7 in [47]. 

Some interesting lines of research appeared in print in 1997-8. First, by keeping the 
cylindrical symmetry, a new diagonal but non-separable family of stiff fluid singularity- 
free solutions was presented in [11]. The family contained the same static limit as the so- 
lution in [27], thereby suggesting that they both form part of a more general class, perhaps 
of non-zero measure, of non-singular cylindrically symmetric stiff fluid models. Second, 
the role of shear in expanding perfect fluid models was analyzed in [6] proving that non- 
singular models with an Abelian spatial G2-symmetry should be spatially inhomogeneous. 
This result was much improved and proven in a more general context in [ 1 ] , showing in par- 
ticular that the symmetry assumption was superfluous. And third, by giving up cylindrical 
symmetry, a family of non-singular (non-perfect) fluid solutions with spherical symmetry 
was presented in [4]. These models depend on one arbitrary function of time. Once again, 
they can avoid the singularity theorems due to the failure of the boundary /initial condition: 
there are no closed trapped surfaces. It was also shown in section 7.8 of [47] that these 
models cannot represent a finite star, since this would require a place where the radial pres- 
sure vanished, which is impossible for appropriate selections of the arbitrary function of 
time. This is a property shared by all models mentioned so far in this subsection. 



1 Recall that all Robertson- Walker geometries are cylindrically symmetric. 
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2.2 Raychaudhuri comes into play; the conjecture 

In December 1995 I attended the International Conference on Gravitation and Cosmology 
(ICGC-95), held in Pune (India), where I had the chance to meet Prof. Raychaudhuri for 
the first time. I was impressed by his personality and accessibility, specially for a man of 
his age and reputation. But more importantly, I was deeply influenced by his remarks in 
brief conversations that we — AKR, Naresh Dadhich and myself — had at that time. If I 
remember well, AKR mentioned averages in these informal conversations, but just by the 
way. This came at a critical time: in a short talk at the workshop on "Classical General 
Relativity", see [5]. I presented the above mentioned combined Robertson- Walker plus 
non-singular general family [46] and its properties. I thought that the paper, which had 
already been accepted, would open the door for realistic models. 

Even though AKR meant spacetime averages, I immediately realized the relevance of his 
idea, especially to discriminate between singular and non-singular cosmological models, 
but using purely spatial averages. As remarked at the end of the previous subsection, all 
known non-singular models were "cosmological" in the sense that they could not describe a 
finite star surrounded by a surface of vanishing pressure. However, it can certainly happen 
that (say) the energy density falls off too quickly at large distances (this certainly occurred 
in all known singularity-free solutions). Thus one may raise the issue whether or not this 
will better describe the actual Universe or rather a weakly-localized object such as a very 
large galaxy. Of course, a good way to distinguish between these two possibilities is to 
use the spatial average of the energy density. Thus, I was inspired by AKR's remarks and 
believed that this was the right answer to the existence of non-singular models such as the 
one in [43]. I incorporated this view to the review [47], see p. 821. 

I met AKR for the second time in Pune again, on the occasion of the 15th International 
Conference on General Relativity and Gravitation (GR15), held in December 1997. Either 
at GR15 or in an informal seminar (I cannot exactly recall), I attended a talk where he 
discussed some of the non-singular models and made some comments about the importance 
of the averages of the physical quantities such as the energy density, the pressure, or the 
expansion of the fluid. In 1998 AKR proved [35] that, under some reasonable assumptions, 
open non-rotating non-singular models must have vanishing spacetime averages of the 
matter and kinematical variables. Later, it was shown [42,48] explicitly that, in the open 
Robertson-Walker models (with an initial singularity), the same spacetime averages vanish 
too. Actually, this holds true for most open spatially homogeneous models as well. Since 
this property is shared by all models, it cannot be used to decide between singularity-free 
and singular spacetimes. 

In my comment, I stressed the following fact. I came to understand, after listening 
to AKR, that pure spatial averages (at a given instant of time) vanish in the known non- 
singular solutions, while they are non-vanishing in open Robertson- Walker models. This, 
together with a well-known singularity theorem for expanding globally hyperbolic models 
(Theorem 1 below), enabled the formulation of the following conjecture [48]: 

In every singularity-free, non-rotating, expanding, globally hyperbolic model 
satisfying the strong energy condition, spatial averages of the matter variables 
vanish. 

This will be made precise and proven in section 4. below. 
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2.3 History of the conjecture 

My comment [48] was only submitted after electronic correspondence with AKR. 2 It 
seems that he was initially skeptical about the use of purely spatial — and not spacetime — 
averages, for he also replied in [36] to our comments. Nevertheless, in an e-mail dated 
on September the 9th, 1998, he mentioned that a letter proving the vanishing of spatial 
averages "following his earlier method" had already been submitted for publication. This 
private announcement was followed, shortly thereafter, by (i) another paper by Dadhich 
and AKR [8] where they proved the existence of oscillatory non-singular models within 
the non-perfect fluid spherically symmetric family of [4] mentioned above, (ii) general 
theorems providing sufficient conditions for the geodesic completeness of general cylin- 
drically symmetric spacetimes [15,12] and (iii) some work [19] showing the relevance that 
the singularity-free solutions might have in the fashionable String Cosmology, see also 
[18] and references therein. 

The letter that AKR mentioned was published in [37]. However, in it global hyperbolic - 
ity was assumed without being mentioned (due to the assumption of the existence of global 
coordinates associated to a hypersurface-orthogonal timelike eigenvector field of the Ricci 
tensor). Moreover, the openness of some local coordinates was taken for granted. Further, 
the statement that the spatial average of the divergence of the acceleration associated to the 
timelike eigenvector field vanishes was not clearly proved. More importantly, the blow-up 
of the kinematical quantities of this eigenvector field was incorrectly related to the blow up 
of some Ricci scalar invariants. 3 All in all, the result in [37], involving spatial averages, 
was not completely proved. 

This led to some interesting works by others [16,17], where the existence of a wide class 
of singularity-free (geodesically complete) cylindrically symmetric stiff fluid cosmologies 
was explicitly demonstrated, and many solutions were actually exhibited. In these papers, 
the family of regular cylindrically symmetric stiff fluids was proven to be very abundant, 
allowing for arbitrary functions. Furthermore, strong support for the conjecture was also 
provided in the second of these papers [17]: the vanishing of the energy density (and 
pressure) of the fluid at spatial infinity on every Cauchy hypersurface was demonstrated 
to be a necessary requirement if the spacetime was to be geodesically complete. This 
was quite encouraging, and constituted the first serious advance towards the proof of the 
conjecture. 

AKR may not have been aware of these important developments and results. He put out 
a preprint [38] to which L. Fernandez-Jambrina [13] found a counter-example. A revised 
version was published in [39] which was also not fully free of the shortcomings, as neatly 
pointed out in [14]. AKR acknowledged these deficiencies in [40], yet however he believed 
that some of his results were still valid. 



2 I thank Dr. Narayan Banerjee for helping with AKR's emails. 

3 Blowing-up of the expansion of the Ricci timelike eigendirection can take place due to many 
different reasons such as a caustic, the failure of hypersurface orthogonality or of local coordinates, 
the failure of the vector field to be a Ricci eigendirection, or a change in the algebraic type of the 
Ricci tensor, among others. There are well-known examples of hypersurface-orthogonal timelike 
congruences showing this behaviour even in flat spacetime. 
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Ref. [40] was AKR's last published paper. In my opinion, after having identified the clue 
to non-singular models (i.e. averages), which led us all to the right track, he tried to prove 
more ambitious and challenging results which were perhaps beyond the techniques he was 
using. This, of course, does not in any way diminishes his fundamental contribution to the 
field of singularities in Cosmology, a subject in which, probably, the most important ideas 
came from his insight and deep intuition. It is in this sense that the spirit of the theorem 
to be proven in section 4. should be credited to him. I can only hope that he would have 
welcomed the new results. 



3. Spatial averages and the Raychaudhuri equation 

3.1 The Raychaudhuri equation 

As is known, the first result predicting singularities under reasonable physical conditions 
was published in 1955 — exactly the year of Einstein's demise — by AKR [33]. In this 
remarkable paper, he presented what is considered to be the first singularity theorem, and 
included a version (the full equation appeared soon after in [24], see also [34]) of the 
equation named after him which is the basis of later developments and of all the singularity 
theorems [21,22,31,47]. The Raychaudhuri equation can be easily derived from the general 
Ricci identity: 

Contracting a with fj, here, then with u v , one gets 

u^V^V^ - u^V^V^ = R pv u p u v 

where R^ v is the Ricci tensor. Reorganizing by parts the first summand on the left-hand 
side, one derives 



i^V^u" + V^Vu" - V„(u"V„u") + R pv u p u v = 



(1) 



which is the Raychaudhuri equation. AKR's important contribution was to understand and 
explicitly show the fundamental physical implications of this simple geometrical relation. 

Let us analyze some of these implications. Observe that in the case that u p defines a 
(affinely parametrized) geodesic vector field, then m^V ' v u p = and the third term van- 
ishes. The second term can then be rewritten by splitting 

V^W;, = S pv + A^ v 

into its symmetric and antisymmetric A pv parts, so that 

W - S^S^ - A^A pv . 

Now the point is to realize two things, (i) If u p is timelike (and normalized) or null, then 
both S^S**" and A^ V A^ V are non-negative, (ii) u p is also proportional to a gradient 



6 



Pramana - J. Phys., Vol. xx, No. x, xxx xxxx 



Singularity theorem 



(therefore defining orthogonal hypersurfaces) if and only if A pv — 0. In summary, for 
hypersurface-orthogonal geodesic time-like or null vector fields u p , one has 

so that the sign of the derivative of the divergence or expansion 8 = V (1 it' i along the 
geodesic congruence is governed by the sign of R pv u p u v '. If the latter is non-negative, 
then the former is non-positive. In particular, if the expansion is negative at some point and 
R pu u p u v > then one can prove, by introducing a scale factor L such that V M (log L) oc 
9 and noting that S^^ = 8, that necessarily the divergence will reach an infinite negative 
value in finite affine parameter (unless all the quantities are zero everywhere). 

If there are physical particles moving along these geodesies, then clearly a physical 
singularity is obtained, since the mean volume decreases and the density of particles will 
be unbounded, see Theorem 5.1 in [47], p. 787. This was the situation treated by AKR for 
the case of irrotational dust. In general, no singularity is predicted, though, and one only 
gets a typical caustic along the flow lines of the congruence defined by u^. This generic 
property is usually called the focusing effect on causal geodesies. For this to take place, of 
course, one needs the condition 

R pv u p u v > (2) 

which is a geometric condition and independent of the particular theory. However, in 
General Relativity, one can relate the Ricci tensor to the energy-momentum tensor T^ v via 
Einstein's field equations (8irG = c = 1) 

R/J.U ~ -jd^vR + Ag^u = Tpv ■ (3) 

Here R is the scalar curvature and A the cosmological constant. Thereby, the condition (2) 
can be rewritten in terms of physical quantities. This is why sometimes (2), when valid 
for all time-like u^, is called the timelike convergence condition, also the strong energy 
condition in the case with A = [21]. One should bear in mind, however, that this is a 
condition on the Ricci tensor (a geometrical object) and therefore will not always hold: see 
the discussion in sect. 6.2 in [47]. 

The focusing effect on causal geodesies predicted by the Raychaudhuri equation was the 
fundamental ingredient needed to derive the powerful singularity theorems. Nevertheless, 
as remarked above, this focusing does not lead to singularities on its own in general. As a 
trivial example, observe that flat spacetime satisfies the condition (2) trivially, but there are 
no singularities: the focusing effect simply leads to focal points or caustics of the geodesic 
congruences. This is why it took some time to understand the necessity of combining the 
focusing effect with the theory of existence of geodesies maximizing the interval (which 
necessarily cannot have focal points or caustics) in order to prove results on geodesic in- 
completeness, which is a sufficient condition in the accepted definition of a singularity. 
With the imaginative and fruitful ideas put forward by Penrose in the 1960s, later followed 
by Hawking, this led to the celebrated singularity theorems, see [31,22,21,47]. 

As a simple but powerful example, which we shall later need to prove the new theorem, 
let us present the following standard singularity theorem (Theorem 9.5.1 in [52], Theorem 
5.2 in [47].) 
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Theorem 1 If there is a Cauchy hypersurface E such that the timelike geodesic congruence 
emanating orthogonal to E has an initial expansion > b > and (2) holds along the 
congruence, then all timelike geodesies are past incomplete. 

The idea of the proof is simple [21,47,52]: since E is a Cauchy hypersurface, the spacetime 
is globally hyperbolic so that one knows that there is a maximal timelike curve from E to 
any point. From standard results any such maximal curve must be a timelike geodesic or- 
thogonal to E without any point focal to E between E and the point. But the Raychaudhuri 
equation (1) implies that these focal points should exist to the past at a proper time less 
than or equal to a fixed value (given by 3/0|s < 3/6.) As every causal curve crosses the 
Cauchy hypersurface E, no timelike geodesic can have length greater than 3/6 to the past. 



3.2 Spatial averages 

Let E be any spacelike hypersurface in the spacetime and let rj s be the canonical volume 
element 3-form on E. The average (f) s of any scalar / on a finite portion S of E is defined 
by 



/ /rfe r 

(f)s = —r = [Vol(S)]- 1 / 

Js 



where Vol(5) is the volume of S C E. The spatial average on the whole E is defined as 
(the limit of) the previous expression when S approaches the entire E 



/ ^ r 

</) E = Jim ^ = [Vol(E)]- 1 J fa (4) 



*7 S 



is 

Obvious properties of these averages are 

1. (linearity) For any S C E, any functions /, g and any constants a, b: 

(af + bg) s = a(f) s + b(g) s 



2. For any S C E, (f) s < (\f\) s . 

3. If S C E is such that its closure is compact (so that its volume is finite Vol(S') < oo), 
then for any / > b > on S, (f) s > b and the equality holds only if / is constant, 
/ = b, almost everywhere on S. In particular, if / > on such an S, then (f) s > 
and the equality holds only if / vanishes almost everywhere on S. 

4. If S C E does not have a finite volume (so that it cannot be of compact closure), then 
for any / > on S, (f) s > and the equality requires necessarily that / — > when 
"approaching the boundary" (i.e., when going to infinity). Conversely, if / > 0, / is 
bounded on S and / is bounded from below by a positive constant at most along a 
set of directions of zero measure, then (f) s = 0. 
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5. If |/| is bounded on S* (|/| < M),then(/ 2 ) s < M(\f\) s . In particular, (|/|) s =0 
implies that (/ 2 ) s = 0. 

6. Similar results hold, of course, for negative and non-positive functions (just use — /.) 



4. The theorem 

Let us consider any spacelike hypersurface S in the spacetime, and let be its unit nor- 
mal vector field (ergo timelike). The projector h^ u = g^ v + u^u v defines the canonical 
first fundamental form of E. A classical result relates the intrinsic Riemannian structure 
of (E, hfj, v ) with the Lorentzian one of the spacetime, and in particular their respective 
curvatures. The main result is the Gauss equation, which reads (e.g. [21,52]) 

Rat3j5hxh^h1h 5 T = Rx^vt + Kx v K^, T - Kx T K^ 

where = K vfJl = h^h^'Vpu^ is the second fundamental form of E, and Rx^vt is 
the intrinsic curvature tensor of (E, h M „). Observe that u^K^ = and u^Rx^vt = 0. 
Contracting all indices here one derives the standard result (e.g. section 10.2 in [52]) 

K 2 = K^K*" + 2R^uV + R-R 

where R is the scalar curvature of S and K = il . Using the Einstein field equations (3) 
this can be rewritten in terms of the energy-momentum tensor as 

K 2 = K^K^ + 2T^u v + 2A-R. (5) 

Recall that T^ v u^u v is the energy density of the matter content relative to the observer w/ 1 , 
and thus it is always non-negative. Note, also, that for any extension of w M outside S as a 
hypersurface-orthogonal unit timelike vector field (still called u M ), their previously defined 
kinematical quantities 6 = V m m m and = V(^u„) are simply 

0\v=K, (S M „ + a( A( M I /))| s = 

where = V ' v u^ is the acceleration vector field of the extended u^. 

Combining this with Theorem 1 one immediately deduces a very strong result concern- 
ing the average of the energy density: 

Proposition 1 Assume that 

1. there is a non-compact Cauchy hypersurface E such that the timelike geodesic con- 
gruence emanating orthogonal to E is expanding and (2) holds along the congru- 
ence, 

2. the spatial scalar curvature is non-positive on average on E: < 0, 

3. the cosmological constant is non-negative A > 0, 

4. the spacetime is past timelike geodesically complete. 
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Then, 

(V«V) S - (K^Kn^ = (#> E = A = . (6) 

Remarks: 

• The first assumption requires spacetime to be globally hyperbolic (so that it is 
causally well-behaved). Also, the timelike convergence condition has to hold along 
the geodesic congruence orthogonal to one of the Cauchy hypersurfaces. Further- 
more, the Universe is assumed to be non-closed (non-compactness) and everywhere 
expanding at a given instant of time — described by the hypersurface E — . All this 
is standard, cf . Theorem 1 . 

• The second assumption demands that the space of the Universe (at the expanding 
instant) be non-positively curved on average. Observe that this still allows for an 
everywhere positively curved S (see the example in the Appendix). This is in accor- 
dance with our present knowledge of the Universe and with all indirect observations 
and measures, e.g. the recent data [53] from WMAP. 

• The third assumption is also in accordance with all theoretical and observational data 
[53]. Observe that the traditional case with A = is included. Notice, however, that 
the timelike convergence condition (2) is also assumed. This may impose very strict 
restrictions on the matter variables if A > 0. 

• The second and third assumptions could be replaced by milder ones such as (i?) < 
2 A, allowing for all signs in both A and The conclusion concerning the 
vanishing of the averaged energy density would be unaltered, as well as the next one 
in (6), but the last two equalities in (6) should be replaced by (-R) s = 2A. 

• Probably, the fourth condition may be relaxed substantially, as one just needs that the 
set of past-incomplete geodesies be not too big. There are some technical difficulties, 
however, to find a precise formulation of the mildest acceptable condition. 

• The conclusion forces the second and third assumptions to hold in the extreme cases 
and, much more importantly, it implies that the energy density of the matter on E 
has a vanishing spatial average. This conclusion was the main goal in this paper. 

Proof: From theorem 1 and the fourth hypothesis it follows that < 9\y, = K cannot 
be bounded from below by a positive constant. Furthermore, the existence of a complete 
maximal timelike curve (which must be a geodesic without focal points) from the Cauchy 
hypersurface E to any point to the past implies that, actually, K can be bounded from 
below away from zero only along a set of directions of zero measure. Point 4 of the list 
of properties for averages implies that (0) s = (A) s = 0, and point 5 in the same list 
provides then ( A 2 ) E = 0. Taking averages on formula (5) and using point 1 in that list 
one arrives at 

(AV^Hs + 2 {T^uV)^ + 2A- (R) E = 0. 

Then, given that all the summands here are non-negative, the result follows. 

This result can be made much stronger by using, once again, the Raychaudhuri equation. 
To that end, we need a Lemma first. 
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Lemma 1 If the energy-momentum tensor satisfies the dominant energy condition and 
(T^u^u")^ = Ofor some unit timelike vector field u^, then al[ the components ofT^ in 
any orthonormal basis {e^} have vanishing average on E 

(T„ v e£e£) E = Ma, = 0, 1, 2, 3. (7) 



Proof: The dominant energy condition implies that [21,49], in any orthonormal basis 
i e a} (where is the timelike leg) 

T„„egeS > |T MV e£e£| , Ma, = 0, 1, 2, 3, 

so that, by taking any orthonormal basis with = eft, points 2 and 6 in the list of prop- 
erties of spatial averages lead to (7) in those bases. As any other orthonormal basis is 
obtained from the selected one by means of a Lorentz transformation — so that the com- 
ponents of T^ v in the new basis are linear combinations, with bounded coefficients, of the 
original ones — the result follows. 

The combination of this lemma with proposition 1 leads to the following result. 

Proposition 2 Under the same assumptions as in proposition 1, ifT^ u satisfies the domi- 
nant energy condition then not only the averages shown in (6) and (7) vanish, but further- 
more 

(«"V M fl - vx> s = («"V„0) E - (iW^) s = 0, (8) 

where v^ 1 is the unit timelike geodesic vector field orthogonal to E, $ = V M t> M its expan- 
sion, is any hypersurface-orthogonal unit timelike vector field orthogonal to E, and 
= u"V ' its acceleration. 

Remarks: 

• The hypersurface-orthogonal vector field may represent the mean motion of the 
matter content of the Universe. Observe that, thereby, an acceleration of the cosmo- 
logical fluid is permitted. This is important, since such an acceleration is related to 
the existence of pressure gradients, and these forces oppose gravitational attraction. 

• The expression u^V^O can be seen as a "time derivative" of the expansion 9: a 
derivative on the transversal direction to E. In particular, v^V^d is the time deriva- 
tive, with respect to proper time, of the expansion for the geodesic congruence or- 
thogonal to E. This proper time derivative of d does have a vanishing average on E. 
Notice, however, that the generic time derivative of the expansion 9 does not have 
a vanishing average in general: this is governed by the average of the divergence of 
the acceleration. 

• Evidently, a^u^ — so that a M is spacelike and tangent to E on E. Furthermore, 
one can write 

VX - {br - uV)V^a v = WVpOv + . 
Letting a represent the spatial vector field a^\^ on E, this implies 

V^a^js = div^a + a- a, 
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where divs stands for the 3 -dimensional divergence within (£, h^ u ) and • is its inter- 
nal positive-definite scalar product; hence a ■ a > and this vanishes only if a = 0. 
Note that the average < div^a > s will vanish for any reasonable behaviour of a M | s , 
because the integral in the numerator leads via Gauss theorem to a boundary (sur- 
face) integral "at infinity", which will always be either finite or with a lower-order 
divergence than vol(E). Therefore, by taking averages of the previous expression 
one deduces 

(VX) S = (a-a)s = W E >0. 

In other words, the first conclusion in (8) can be rewritten in a more interesting way 

as (u^Vfj.O — a v a v )^ = 0, or equivalently 

(u"V„0) 2 = (a,a") E > 0. 

Observe that, if these averages do not vanish, this implies in particular that there 
must be regions on E where u^V^O is positive. 

Proof: From proposition 1 it follows that necessarily A = 0, hence in any orthonormal 
basis {e^} the equations (3) imply 

Rn^ e a e = T^e^ep — -T p p r) a/ 3, 

where {r} a fj) =diag (—1, 1, 1, 1). By taking averages on £ here and using lemma 1 one 
deduces 

(iVe£e£) E = Va, = 0, 1,2, 3 

which are the last expressions in (8). 

The Raychaudhuri equations (1) for and are respectively 

u» \7 U 9 + V^Vu" - VX + R^u v = 0, (9) 
v v V v ti + V^VV + R^v v = 0. (10) 

Obviously u M |s = w^ls and it is elementary to check that 9\s = #|e and 

V^Vu^s = V^WIs = K^K^. 
From (9-10) one thus gets 

{u v VJ - V„o")Ie = w"V v i?| s = -K^K^ - R^u v \^ < 

so that u v V v 6 — V M a M is everywhere non-positive on S. Taking averages on S, using the 
second in (6) and the previous result on the averages of the Ricci tensor components one 
finally gets 

(u"V„fl - V M a">ls = (u"V^)| E - 0. 

That ends the proof. 

The main theorem in this paper is now an immediate corollary of the previous proposi- 
tions. 
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Theorem 2 Assume that 

1. there is a non-compact Cauchy hypersurface E such that the timelike geodesic con- 
gruence emanating orthogonal to E is expanding and the timelike convergence con- 
dition (2) holds along the congruence, 

2. the spatial scalar curvature is non-positive on average on E: (i?) E < 0, 

3. the cosmological constant is non-negative A > 0, 

4. the energy-momentum tensor satisfies the dominant energy condition. 
If any single one of the following spatial averages 

A, <0> E , <0) E , <fl 2 ) s , (tf 2 ) E , {K I1V K I1V ) T . , (i?) E , (v"V„0) E , 
W-V/) E , (T^e£e£) E , (iV<e£) s Va, = 0,1, 2, 3 

(foes nof vanish, then the spacetime is past timelike geodesically incomplete. 

As before, notice that the first average in the second row can be replaced by 
(u^V^O — a^a M ) E . Therefore, another sufficient condition is that u^V^O is non-positive 
everywhere but its average be non-zero. 

5. Conclusions 

Observe the following implication of the theorem. Under the stated hypotheses, a non- 
vanishing average of any component of the energy-momentum tensor (or of the Ricci ten- 
sor) — such as the energy density, the different pressures, the heat flux, etc. — leads to the 
existence of past singularities. It is quite remarkable that one does not need to assume any 
specific type of matter content (such as a perfect fluid, scalar field, ...); only the physi- 
cally compelling and well established dominant energy condition is required. The theorem 
is valid for "open" models, as the Cauchy hypersurface is required to be non-compact. 
This is, however, no real restriction, because for closed models there are stronger results 
[21,22,31,52,47]. As a matter of fact, closed expanding non-singular models require the 
violation of the strong energy condition [47]. There is no hope of physically acceptable 
closed non-singular models satisfying the timelike convergence condition (2) (without this 
condition, there are some examples [29,47]). 

Let me stress that the conclusion in theorem 2 is quite strong: it tells us that the incom- 
pleteness is to the past. Besides, I believe that one can in fact prove a stronger theorem 
such that the geodesic timelike incompleteness is universal to the past. 

The main implication of the above results is this: a clear, decisive, difference between 
singular and regular (globally hyperbolic) expanding cosmological models is that the latter 
must have a vanishing spatial average of the matter variables. Somehow, one could then say 
that the regular models are not cosmological, if we believe that the Universe is described 
by a more or less not too inhomogenous distribution of matter. This is, on the whole, a 
very satisfactory result. 
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Appendix: the singularity-free model of [43] 

In cylindrical coordinates {t, p, tp, z] the line-element reads 

ds 2 = cosh 4 (ai) cosh 2 (3a p)(-dt 2 + dp 2 ) + 
^2 cosh 4 (at) coslT 2/3 (3ap) sinh 2 (3ap)aV 2 + cosh~ 2 (ai) cosh~ 2/3 (3ap)dz 2 , (11) 

where a > is a constant. This is a cylindrically symmetric (the axis is defined by p — > 0) 
solution of the Einstein field equations (3) (with A = for simplicity) for an energy- 
momentum tensor describing a perfect fluid: T^ v = gu^Uy + p (g^ u + u^u v ). Here g is 
the energy density of the fluid given by 

q = 15a 2 cosh -4 (at) cosh _4 (3ap) , 

and 

Up = (— cosh 2 (ai) cosh(3a/9), 0, 0, 0) 

defines the unit velocity vector field of the fluid. Observe that is not geodesic (except 
at the axis), the acceleration field being 

a M = (0,3atanh(3ap),0,0) . 

The fluid has a realistic barotropic equation of state relating its isotropic pressure p to g 
by: 

1 

This is the canonical equation of state for radiation-dominated matter and is usually as- 
sumed to hold at early stages of the Universe. Note that the energy density and the pressure 
are regular everywhere, and one can in fact prove that the space-time ( 1 1 ) is completely free 
of singularities and geodesically complete [2]. For complete discussions on this spacetime 
see [2] and sect. 7.6 in [47]. One can nevertheless see that the focusing effect on geodesies 
takes place fully in this space-time. This does not lead to any problem with the existence of 
maximal geodesies between any pair of chronologically related points, see the discussion 
in [47], pp. 829-830. 

Space-time (11) satisfies the strongest causality condition: it is globally hyperbolic, 
any t =const. slice is a Cauchy hypersurface. All typical energy conditions, such as the 
dominant or the (strictly) strong ones (implying in particular (2) with the strict inequality) 
also hold everywhere. The fluid expansion is given by 
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fl = VX = 3a 3 Sinh(at) • (12) 

cosh (at) cosh(3ap) 

Thus this universe is contracting for half of its history (t < 0) and expanding for the second 
half (t > 0), having a rebound at t = which is driven by the spatial gradient of pressure, 
or equivalently, by the acceleration a^. Observe that the entire universe is expanding (that 
is, 8 > 0) everywhere if t > 0. Note that this is one of the assumptions in propositions 
1, 2 and theorem 2. It is however obvious that, for any Cauchy hypersurface St given by 
t = T =constant, the average < 8 >s T = 0. As one can check for the explicit expression 
(12), 9 is strictly positive everywhere but not bounded from below by a positive constant 
because lim^oo 8 = 0. Observe, however, that for finite p one has lim z ^oo 8 > and 
finite. 

Similarly, one can check that the scalar curvature of each St is given by 

R = 30a 2 cosh" 4 (aT) coshT 4 (3ap) > 

which is positive everywhere. However < R >s T ~ 0, and analogously < g >s T =< 
V >e t = 0. Observe also that 



sinh 2 (3op) 



2 



a M a M = 9a 

cosh (at) cosh (3ap) 

and thus < a M a M >s T = 0. This implies that in this case < u^V ^9 >£ T = 0. The sign of 
m m V ^9 — V^a^ 1 is negative everywhere, as can be easily checked: 



cosh (at) cosh (3ap) cosh (at) cosh (3ap) 

All this is in agreement with, and illustrates, the main theorem 2, propositions 1 and 2, and 
their corresponding remarks. 

This simple model shows that there exist well-founded, well-behaved classical mod- 
els which expand everywhere, satisfying all energy and causality conditions, and are 
singularity-free. However, as we have just seen, the model is somehow not "cosmolog- 
ical" to the extent that the above mentioned spatial averages vanish. 
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